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We take both the bilinear and trilinear R-parity violating couplings in supersymmetric models 
as a source of neutrino masses and mixings. Using the solar and atmospheric data and the 
Chooz constraint we determine the allowed ranges of those couplings. We also estimate the 
effective mass for neutrinoless double beta decay in this scenario. 



1 R-parity violation and its parametrization 

In supersymmetric models, lepton and baryon numbers {L and B, respectively) are not auto- 
matically conserved, unlike in the standard model. In terms of L and B, one defines R-parity 
as R p = (— i) 3 - B + L + 2S ' ; where S is the spin of a particle. Stringent phenomenological limits on 
R-parity violating ($ p ) interaction^^ can be found in Ref.^. B violation has got nothing to 
do with neutrino mass generation, and we assume that such interactions are absent. Neutrino 
Majorana mass generation requires two units of L violation. For this purpose, we allow both 
bilinear (/ij) and trilinear (A, A') interactions in the superpotential as well as the bilinear soft 
terms (Bi), given by, 

W = n J H u Lj + ±\ JKt LjL K El + >t Jm LjQ p Dl + h™H u Q v U c q , (1) 
where the vector Lj = {H^Lj) with J : 4..1. The soft supersymmetry-breaking potential is 

V soft = ^HtH u + h^[m 2 L }j K L K + B J H u Lj 

+A u p s H u Q p U c s + A j p s LjQ p D c s + ±A JKl LjL K Ef + h.c. (2) 

It should be noted that field redefinitions of the Ha, Li fields correspond to basis changes 
in Lj space and consequently the Lagrangian parameters will be altered. Hence we use the 
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basis-independent parameters constructed in™" and write the neutrino mass matrix in terms of 
such parameters S^Sg, , 5y q , which in the basis in which the sneutrino vacuum expectation 
values are zero correspond to the Lagrangian parameters n l B l /\B\, Xijk, AL g , respectively. 

This talk is main ly ba sed on the paper^. For previous studies of Ifip effects on neutrino 
masses, see also Refs.™^. 



2 Generation of neutrino masses by R-parity violation 

To understand the effects of the bilinear and trilinear terms on neutrino mass matrices, let us 
consider them one by one. First, switch on only the bilinear 5^ terms which appear in the 
superpotential. Such terms generate neutrino masses at the tree level which are proportional to 
5^5^. This constitutes a rank 1 mass matrix which leads to only one nonzero eigenvalue. This 
is not enough since we know that the solar and atmospheric neutrino data require at least two 
nonzero eigenvalues. 

Then we turn on the bilinear 5 B terms which appear in the scalar potential. They contribute 
to neutrino mass at the loop level via the Grossman-Haber diagram in which there are gauge 
couplings at the neutrino vertices while there are two types of Ijip interactions giving rise to the 
AL = 2 Majorana mass via the diagram of Fig. ^ The first type has Ijip couplings located at 
positions III + IV (slepton-Higgs mixing) with contributions proportional to 5b$b- in the second 
type, the Ijip interactions are located at positions V + IV (neutrino-neutralino and slepton-Higgs 
mixing) with contributions proportional to 5^5b- Now, the 8 % „ and 5 B terms together give rise 
to two nonzer o ma sses . while one eigenvalue still remains zero. As far as data is concerned, this 
scenario works^^. 

Then we turn on the trilinear interaction for L violation, namely, the 5 l ^ k (or 5y q ) couplings. 
We then have loops involving the trilinear Ijip couplings A or A' at the neutrino vertices I and 
II in Fig. 1 (with lepton/slepton or quark/squark as propagators). They give contributions 
proportional to b~\o~\ (or 6\'S\i). One also has the diagram of Fig. [21 where two units of L 
violation come from positions V (neutrino-neutralino mixing) and II (A or A' vertex). Their 
contribution to the neutrino mass is proportional to 5^5\ ( or 8u8\')- Now, with bilinear and 
trilinear terms together, all the neutrinos become massive. 
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Figure 1: The usual loops (Iji v at I + II) and the Grossman-Haber loops (_$ p at III + IV or V + IV) contributing 

to the neutrino mass. 

For the sake of simplicity, we set all unknown sparticle masses equal to M$ = 100 GeV. We 
then obtain a neutrino mass matrix of the form 

[m v }ij = M s 
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Figure 2: Loops with a gauge and a trilinear Yukawa coupling. The interactions are located at V and II. 
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We have included for completeness the contributions arising from the 5\> terms which we set to 
zero in our numerical analysis. The simplest case to consider with a common 5^ = <5 M , 5 B = 5b 
and 5™ k = 5\ does not work as it cannot accomodate two large mixing angles for B\i and #23- 

In our numerical analysis we take : 5^, S l B , 5™ k = 5\, 5™ k = 0, for % = 1, 2, 3, i.e., seven 
independent parameters. Note that a common 5\, in addition to 5^ and 5 l B terms, is enough 
to make all the neutrinos massive, and allows us to study how the presence of the trilinear 
interaction alters the allowed range of the bilinear parameters. Thus, the neutrino mass matrix 
elements will be given by, 
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where we have employed the hierarchy of the charged fermion masses to keep only the dominant 
terms, i.e. those which involve m T . 



3 Experimental data on neutrino masses and mixings 



We write the PMNS matrix, which is the rotation matrix from neutrino flavour (/) to mass (i) 
eigenstates, as 



V, 



fi 



C12C13 C13S12 S13 

-C23S12 - C12S13S23 C12C23 - S12S13S23 C13S23 
S23S12 - C12C23S13 -Cl2«23 - C23S12S13 C13C23 



(6) 



where Cjj = cos 9ij and = sin 9ij . We have neglected the CP phases which are not relevant for 
our purpose of extracting allowed ranges of new physics from oscillation analysis. We take the 
solar anomaly to be a consequence of v e -Vn oscillation, and the relevant mass squared difference 
is Am 2 2 = Am 2 olar . We consider the atmospheric oscillation to be between and u T , and 
the relevant mass squared difference is Am| 3 ~ Am| 3 = Am 2 tm . After the announcement 
of the SNO data, the MSW-LMA oscillation is the most favoured solution with Am^ lar = 
(2.5 — 19.0) x 10~ 5 eV 2 and sin 2 2#i 2 = 0.61 — 0.95. The SuperK atmospheric neutrino data 
suggest Am 2 tm = (2 - 5) x 10~ 3 eV 2 with sin 2 26> 2 3 = 0.88 - 1.0. The Chooz and Palo Verde 
long baseline reactor experiments constrain sin 2 #13 ^ 0.04. Tritium /3-decay requires that the 
absolute mass is m Ve 2.2 eV. For the references of all the experimental data, see" 

4 Results and conclusions 



We have performed a general scan of parameter space made up by the seven parameters 
(three 5* three 5 l B , one 5\) that appear in the mass matrix. We allow the tree-level 
contributions to either dominate over the loop corrections, to be on the same order as 
these, or to be much smaller than the loop terms. The fitted values of the couplings, 
satisfying all the data mentioned in the previous section, are given in table I. 



In Fig. |3] we have presented the allowed region in the \6b\ = v/Eife) 2 versus \5 t 



J2i($li) 2 plane for the combined fit. We show our results for both 5\ ^ and 5\ = 0. It 
can be seen that the allowed region increases when we admit non-zero values of 5\. This 
happens due to the presence of the 5^5\ terms in the mass matrix (originating from Fig.[2J) 
which can take either sign thus accomodating a larger region of the parameter space. This 
figure should be compared with Fig. 5 of Ref.!* 
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The resulting fit strongly prefers a hierarchical mass pattern in our scenario, although a 
distinction between the inverted and the normal hierarchy is not possible. 

The maximum value of m e s we have predicted (see Table 1) can hopefully be tested in the 
next generation of neutrinoless dou ble beta decay experiments. The range again implies 
that the spectrum is not degenerate^! 

Allowing a non- vanishing 5y will not qualitatively change the pattern of our fit. 

The analysis was done before the first announcement of the KamLAND results. The 
pattern of the fit will not be qualitatively altered if we include the KamLAND data. 
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Couplings 


Min 


Max 




-2.0 x 10~ 4 


2.0 x 10~ 4 


si 


-6.8 x 10-' 


6.8 x lO" 7 


si 


-8.4 x 10-' 


8.4 x 10-' 


si 


-8.4 x 10~ v 


8.4 x lO" 7 




-2.7 x 10~ 5 


2.7 x 10~ 5 


s% 


-3.0 x 10~ 5 


3.0 x 10~ 5 


s% 


-3.0 x 10~ 5 


3.0 x 10~ 5 


m cfr (eV) 


1.9 x 10~ a 


6.7 x 10~ 2 




4.9 x 10~ 2 


0.2 



Table 1: Allowed range of the couplings satisfying MSW-LMA, SuperK and Chooz simultaneously (with cos (3 — 1). 
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Figure 3: We present the allowed region in the 8b = \Sb\ versus <5 M = |5 M | plane. The circles are solutions for 

<5a 7^ 0, the crosses are for 5\ = 0. 



